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ABSTRACT 

We present IRAM Plateau de Burc Interferometer observations of the 12 CO (3-2) 
' emission from two far-infrared luminous QSOs at z ~2.5 selected from the Herschel- 

j_j , ATLAS survey. These far-infrared bright QSOs were selected to have supermassive 

Ci i black holes (SMBHs) with masses similar to those thought to reside in sub-millimetre 

galaxies (SMGs) at z ~ 2.5; making them ideal candidates as systems in transition 
from an ultraluminous infrared galaxy phase to a sub-mm faint, unobscurcd, QSO. We 
detect 12 CO (3-2) emission from both QSOs and we compare their baryonic, dynamical 
and SMBH masses to those of SMGs at the same epoch. We find that these far-infrared 
bright QSOs have similar dynamical but lower gas masses than SMGs. In particular 
we find that far-infrared bright QSOs have ~50±23% less warm/dense gas than 
SMGs, which combined with previous results showing the QSOs lack the extended, 
cool reservoir of gas seen in SMGs, suggests that they are at a different evolutionary 
stage. This is consistent with the hypothesis that far-infrared bright QSOs represent a 
short (<~ 1 Myr) but ubiquitous phase in the transformation of dust obscured, gas-rich, 
starburst-dominatcd SMGs into unobscured, gas-poor, QSOs. 

Key words: galaxies: formation, — galaxies: evolution — quasars: emission lines — 
quasars: individual — J0908-0034, J0911+0027 
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1 INTRODUCTION 

An evolutionary link between local ultraluminous infrared 
galaxies (ULIRGs, with far-infrared [FIR] luminosities of 
> lO^ La) and QSOfl was first suggested by 
I Sanders et al.l (| 1988h and considerable effort has been ex- 
pended to test this connection in the lo cal Universe (e.g. 
iTacconi et ail |2002| ; rVeilleux et all 120091 ). This hypothesis 
has been strengthened by the discovery of a relation be- 
tween the mass of supermassive bla ck holes (SMBHs) and 
the mass of their host spheroids (e.g. Magprrian et al.ll 19981 ; 
iGebhardt et al.ll200d ; lFerrarese fc MerrittboOOT ). which sug- 
gests a physical connection between the growth of spheroids 
and their SMBHs. As shown by theoretical simulations, such 
a link can be formed through the suppression of star forma- 
tion in a host galaxy by winds and o utflows from an active 
galactic nucleus (AG N) at its centre (|Di Matteo et al1l2005l ; 
iHopkins et ai1l2005l ). 

Both the ULIRG and QSO populations evolve very 
rapidly with redshift, and although subject to selec- 
tion biases, intriguingly both populations appear t o reach 
a br o ad peak in their ac ti vity at 2 ~ 2.5 (IShaverl 
1 19881 ; IChapman et~aH 120051 ; IWardlow etHI 120111 ). The 
high-redshift ULIRGs, which are typically bright at sub- 
millimetre wavelengths, and hence are called sub-millimetre 
galaxies (SMGs), are thought to represent the formation 
of massive stellar systems through an intense burst of 
star formation in the early Universe, triggered by merg- 
ers of gas-rich progenitors dFraver et ahl 19981; [Blain et alj 
2002l:ISwinbank et al.l20M l2010al ; lTacconi et al.l2006l . l200Sl ; 



Engel et al.ll2010l ) 



The similarity in the redshift distribution of the SMGs 
and QSOs may be indicating that the evolutionary link 
between these populations, which has been postulated lo- 
cally, also holds at z ~ 2, when both populations were 
much more numerous, and significant in terms of stellar 
mass assembly. Indeed recent work on the clustering of 
z ~ 2 QSOs and SMGs has shown that both populations 
reside in parent ha los of a similar mass, Mhaio ~ 10 13 M© 
(|Hickox et al.ll2012l ). adding another circumstantial connec- 
tion between them. Moreover, this characteristic Mhaio is 
similar to the mass at which galaxy pop ulations transition 
from star forming to passive systems (e.g. lBrown et af 1l200Sl ; 
ICoil et al.lboOSl ; Hartley et alj|2010l ) which may be indirect 
evidence for the influence of QSO-feedback on star f orma- 
tion in galaxies l|Granato et al. l200ll ; lLapi et al.|[20 06). One 
other piece of circumstantial evidence for a link between 
SMGs and QSOs comes from estimates of the masses of 
SMBHs in SMGs. Due to the dusty nature of SMGs these 
studies are inherently challenging, yielding results with more 
scatter than seen for optical quasars. However, they suggest 
that the SMBH masses are ~ 10 8 Mq , significantly lower 
than seen in comparably m assive galaxies at the present- 
day (| Alexander et al]|2008l ). This indicates the need for a 
subsequent phase of SMBH growth, which could be associ- 
ated with a QSO phase. 

One puzzling issue about the proposed evolutionary 
connection between high-redshift QSOs and SMGs is only a 
small fraction of sub-millimeter detected galaxies (Ssso ^ 

1 QSO, quasi-stellar object, defined as having Mb < —22 and 
one or more broad emission lines with a width > lOOOkms -1 
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Figure 1. The 2S LAQ UV rest frame spectra for both QSOs 
in our sample from [Croom et al . (2009), with the 2 QZ compos- 
ite QSO spectrum overlaid (red; [Croom et al]|2002 i). From left- 
to-right, a dashed line indicates the rest wavelength of Lyct, 
Nv, Siiv, Civ and Cm] emission. Both QSOs, but especially 
J0911+0027, have relatively weak Lyc? 1215 emission when com- 
pared to Civ 1549, suggesting the presence of sig nificant quan- 
tities of neutral gas in their vicinity (see also lOmont et alj 
Il996l) . For each QSO we estimate the SMBH mass from the 
FW HM of the Civ 1549 line and the rest frame 1350A luminos- 
ity flVestereaard fc Peterson! |2006| ) . The FWHM of the Civ line 
is derived from the best fitting Gaussian and continuum model 
(shown in blue), which provides an adequate fit to the emission 
line in both cases: x? = 1-4 for J0908-0034 and x? = 1-5 for 
J 091 1+0027. The rest wavelength scale of the plots is based on 
the systemic redshifts derived from 12 CO (3—2), see section §[3] 



5mJy), are identified as optically-luminous, Mb < —22, 
z ~ 1 -3 QSOs (~ 2%. IChapman et al.ll2005l ; IWardlow et all 
l201ll ). indicating that, if they are related, then the QSO and 
ULIRG phases do not overlap significantly. Taken together, 
these various results provide support for an evolutionary 
link between high-redshift ULIRGs (SMGs) and QSOs, 
with a fast transition (~ lMyr) from the star-formation- 
dom inated SMG-pha se to the AGN-dominated QSO-phase 
(e.g. IPaee et al.ll2012h . where rapid SMBH growth can then 
subsequently account for the present-day relation between 
spheroid and SMBH masses. In this scenario, when a QSO 
is detected in the far-infrared/sub-millimetre it will be in 
the transition phase from an SMG to an unobscured QSO, 
making its physical properties (e.g. gas mass, dynamics, etc.) 
a powerful probe of the proposed evoluti onary cycle (e.g. 
IPage et al.ll2004 l201ll ; IStevens et al]|2005f ). 

In this paper we present a study with the IRAM Plateau 
de Bure interferometer (PdBI) of the cold molecular gas 
in far-infrared-selected "transiti on" candidate QS Os from 
the Hersche^ H- ATLAS survey (|Eales et al.ll2010l ). We de- 
tect 12 CO (3-2) in both far-infrared bright QSOs studied, 
and compare the gas and kinematic properties of these 
to other 12 CO-detected QSOs and SMGs to relate their 
evolu tion. Throughout w e adopt cosmological parameters 
from ISpergel et ail (|2003l ) of: U m = 0.27, Q. A = 0.73 and 
H = 71kms- 1 Mp C - 1 . 



2 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA 
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Figure 2. Upper. The velocity-integrated 12 CO (3-2) emission for the two QSOs (summed over the regions marked in the spectra below). 
The contours represent signal-to- noise, and are spaced at ±3, 4, 5, ... X a (dashed contours are neg ative). Significant 12 CO(3-2) emission 
is detected from both QSOs and is well-centered on the optical positions of the QSOs, indicated by a cross (offset by 1-2", which is within 
the primary beam). Analysis of the channel maps shows both sources are unresolved with no evidence of velocity gradients, although 
the modest spatial resolution of our maps means that this is not a strong constraint. Lower. Spectra showing the redshifted 12 CO (3—2) 
emission lines from both QSOs. The spectra were binned to a resolution of 10 MHz and 20 MHz respectively, and we overlay the best 
fitting Gaussian plus continuum model (the continuum emission in both QSOs is negligible). The redshift derived from 12 CO (3-2) as 
well as the line flux and line widths derived from the Gaussian fits are reported in Table 2. We note that the FWHM of the J0908— 0034 
emission line is relatively narrow, at 125±25kms — 1 , compared to 530±100kms _1 for J0911+0027. In each spectrum the predicted 
postion of the 12 CO (3-2) emission line, from the 2SLAQ UV spectrum, is indicated by an arrow. The two dot-dashed lines correspond 
to ±FWHM, and show the region which was summed to produce the maps above. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Sample Selection 

Our aim is to study the molecular gas and dynamical prop- 
erties of far-infrared bright QSOs to test their evolutionary 
link to SMGs, for which larger sampl es are now available 
(|Greve et al.l2005l : lBothwell et al.ll2012l ). We focus our study 
on potential transition (far-infrared bright) QSOs at the era 
when the activity in both the QSO and SMG populations 
peaked, z ~ 2.5. Previous 12 CO stud ies of unlensed far - 
infrared or sub-millimetre-bright QSOs (jCoppin et al.|[2008l ) 
have usually studied the most luminous QSOs, with SMBH 
masses 1-2 orders of magnitude larger than those seen in 
typical SMGs and correspondingly low space densities. Both 
of these factors make it hard to relate the properties of these 
extr eme QSOs to ty pical SMGs (studies of lensed QSO sam- 
ples, iRiecherslliHll . suffer from the difficulty of determining 
the true space densities of the sources being studied). The 



limitations of these previous studies arose because of the 
lack of l arge samples o f far-in frared /sub-millimetre detected 
QSOs; ICoppin et ail (|2008l ) selected a sample of QSOs 
which were later found to be bright in the sub-mm. But with 
the advent of wide-field far-infrared / s ub-millimetre survey s 
with the Herschel Space Observatory (jPilbratt et a.1.1 12010T ) 
we can now select samples of QSOs with redshifts, space 
densities, far-infrared luminosities and most critically SMBH 
masses which are well-matched to those of the hypothesised 
descendents of typical SMGs. 

For our analysis we have there fore used the s ample 
of far-infrared observed QSOs from iBonfield et ail (|201ll ) 
whic h are derived from the survey of the 9 h H-ATLAS 
field (|Eales et al.ll201fj| ). lBonfield et all (|201ll ) extracted far- 
infrared fluxes from the SPIRE 250, 35 and 500-/im maps 
i|Grifnn et al.ll201fj| ; iPascale et aill201ll) at the optical p osi- 
tion of 3 72 QSOs from eithe r SPSS (|Schneider et al.l2010h or 
2SLAQ (|Croom et al.ll2009h surveys across the ~ 16 degree 



© 0000 RAS, MNRAS 000, 000-000 



4 Simpson et al. 



Table 1: Observed Properties 



QSO 


"CO 


<5co 


S 250/jm a 


S 350/im a 


S 500/jm a 


S 3mm 




(J2000) 


(J2000) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


J0908-0034 


09 08 47.18 


-00 3416.6 


9.0±6.6 


26.2±8.1 


15.0±9.0 


< 1.4 


J0911+0027 


091148.30 


+00 2718.4 


25.3 ±5.8 


14.7±6.6 


7.7 ±7.8 


< 1.0 



a J0908-0034 flux density values are after deblending sources. Fluxes extracted at the optical position of the QSO, before deblending, 
are 15.6 ± 6.4 mjy, 39.5 ±7.1 mjy and 41.0 ±8.5 mjy 



field. A statistical background subtraction was performed on 
the measured flu xes for each QSO to a ccount for confusion 
in the maps, see iBonfield et all (|201ll ). Of the 372 QSOs, 
~ 20 % are consistent with having Lfir > 10 12 Lo, mak- 
ing them candidates for "transition" QSOs (e.g. between a 
far-infrared bright SMG and far-infrared faint unobscured 
QSO phase). Restricting the redshift range to z = 2-3 
(to roughly match the redshift peak of SMGs, Chapman 
et al. 2005), gives 29 QSOs. We then estimated SMBH 
masses for these QSOs from the extrapolated rest frame 
5100-A luminosities (e.g. Wandel et al. 1999) to select those 
with S MBH masses of M B h < 3 x 10 s M , comparable to 
SMGs (I Alexander et all [2008). These selection criteria re- 
sulted in nine targets and from this list we selected six 
QSOs for a pilot project to determine their gas masses 
and dynamics through the detection of 12 CO (3-2) emis- 
sion with PdBI. We estimated systemic redshifts for these 
QSOs using the measured wavelengths from Gaussian fits 
to the Siiv 1397, Civ 1549 and Cm] 1909 emission lines in 
their spectra (Fig. [T]). We then adjust our estimates by the 
weighted mean velocity offsets of t hese lines from 12 CO (3- 
2), for the 12 CO-detected QSOs in lCoppin et~aH (|200Sf ) (we 
find Siiv 1397 provides the best estimate of the systemic 
redshi fts). We note that, as previously seen bv lOmont et alj 
(| 19961 ). roughly ~ 80% of these far-infrared bright QSOs 
show very weak Lya emission. There may be a relation be- 
tween far-infrared luminosity and weak Lya, but its physical 
origin is still unknown. 

This paper presents the PdBI observations of the first 
two QSOs from this sample: J0908-0034 (RA: 09 08 47.18, 
Dec: -00 3417.9, J2000; z vv = 2.5073) and J0911±0027 
(RA: 091148.38, Dec: ±00 2719.7, J2000, z vv = 2.3697). 



2.2 PdBI Observations 

We used the six-element IRAM PdBI in compact (D) con- 
figuration to search for redshifted 12 CO (3-2) emission from 
the QSOs. We tuned the correlator WideX to the frequency 
of redshifted 12 CO (3-2) from the estimated redshifts de- 
rived from their UV emission lines (Fig. [TJ: 98.593GHz 
and 102.619GHz, for J0908-0034 and J0911±0027 re- 
spectively. The advantage of WideX is that its 3.6-GHz 
(dual polarisation) spectral coverage, at a fixed channel 
spacing of 1.95MHz, corresponds to a velocity range of 
~ 10,000 kms -1 , sufficient to identify the 12 CO(3-2) line 
even if the UV lines are significantly offset from the sys- 
temic redshift. We obtained a total on-source observing time 
of 6hr per source between 2010 April 26 and 2010 May 4. 
The overall flux density scale was calibrated to 3C 273, with 
observations of J0906±015 for phase and amplitude cali- 



bration. Receiver bandpass calibration was performed on 
0923±392. The data were calibrated in the GILDAS software 
package and a naturally weighted data cube produced for 
each QSO. In this configuration the synthesised beam for 
natural weighting is 6.3" x 5.0" at a position angle (PA) of 
166° for J0908-0034, and 7.3" x 5.6" at a PA of 19° for 
J0911±0027. 

2.3 Far-infrared Luminosities and SMBH Masses 

For our analysis we require more accurate estimates of the 
far-infrared luminosities and SMBH masses of our QSOs. To 
estimate far-infrared luminosities of the QSOs in our sam- 
ple, we exploit the Herschel SPIRE 250, 350 and 500-pm 
imaging. For each QSO, we deblend multiple sources in the 
vicinity, using the 250/^m map to identify nearby sources. 
For J0911±0027 there is a bright 250/im source centered at 
the optical position of the QSO which is clearly isolated, but 
J0908-0034 has significant ( 40%-60% at 250-500/xm) con- 
tamination from a galaxy ~30" away. We report the final 
deblended photometry for each source in Table 1. To derive 
the far-infrared luminosity for each QSO, we fit a modi- 
fied black body spectra to the SPIRE photometry, adopt- 
ing a dust emissivity, j3 — 1.6 and fixing the temperature 
of the modified black body at T<j = 40 K, (the averag e 
dust temperature of z ~ 2-3 QSOs, iBeelen et al.l 120061 ) 
and calculate the far-infrared luminosity (Lpir) by integrat- 
ing the best-fit SED between (rest-frame) 8-1000 ^m. We 
note that if we instead allow the characteristic dust tem- 
perature to vary in the SED fit, we see at most a ~45% 
change in Lpir, consistent within our error estimates. The 
modified black body fits confirm the Lpir for both QSOs is 
~ (3 - 4 ± 1) x 10 12 L (Table 2). The QSOs are thus sig- 
nificantly fainter in the FIR than the QSO sample of Cop- 
pin et al. (2008), Lpir = (7.5 ± 1.5) x 10 12 L Q , and also 
at the fainter end of the SMGs in Bothwell et al. (2012), 
Lfir = (4.8 ± 0.6) x 10 12 L© . The implied star-formation 
ratef0 (SFR) are 500-700 M Q yr" 1 . 

We calculated SMBH masses for both QSOs from the 
FWHM of their Civ 1549 emission lines and their rest frame 
1350- A continuum luminosit i es (L 1350), following Eq. 7 
from IVestergaard fc Peterson! (|2006l ). This relation is cal- 
ibrated to the results of reverberation mapping, however 
the geometry of the broad line region (producing the emis- 
sion lines) is poorly con strained and may bias virial mass 
estim ates to low values (|jarvis fc McLurdkood : I Fine et al.l 
l201lf ). We derive the FWHM from a Gaussian fit to the 

3 SFR (Meyr" 1 ) = 1.7 x 10- 10 L PIR (L ) iKennicutil l ll998t) . 
following a Salpeter IMF over a mass range 0.1-100Mq 
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Table 2: Physical Properties 



Source 


zco 


ScoAi/ 
(Jy kins" 1 ) 


FWHMco 
(kms -1 ) 


T / 

^CO (3-2) 
(10 10 Kkms" 1 pc 2 ) 


Lfir 
(10 12 L Q ) 


M gas 
(10 10 M o ) 


M BH 
(10 s M ) 


J0908-0034 
J091 1+0027 


2.5008 ±0.0002 
2.3723 ±0.0005 


0.23 ±0.04 
0.62 ±0.10 


125 ±25 
530 ± 100 


0.77 ±0.15 
1.87 ±0.33 


3.1±0.8 
3.5±0.8 


0.77±0.15 
1.9±0.3 


1 6 +0 ' 3 

2 5+ 1 ' 4 
z -°-0.9 



ClV emission line in the 2SLAQ spectra (Fig. [T]), measuring 
4500±340kms" 1 and 5900±i|go kms" 1 for J0908-0034 and 
J0911±0027 respectively, and determine L1350 from the g- 
band SDSS magnitudes of g = 21.41 ± 0.05 for J0908-0034 
and g = 21.64 ± 0.05 for J0911±0027. As we do not see 
strong evidence of reddening in the 2SLAQ spectrum of ei- 
ther QSO we do not consider the effects of dust extinction 
on our estimates of L1350 (Fig- [I]). From these we estimate 
SMBH masses of ~ 2 x 10 8 M for both QSOfQ (Table 2). We 
caution that estimates from the Civ emission line can over 
predict the SMBH mas s by a factor of 2-5, when compared 
to estimates from Ha (|Ho et al.ll2012T ). In the absence of a 
more robust mass estimator, we simply acknowledge that 
the SMBH masses of our sources may be over estimated. 



3 ANALYSIS & RESULTS 

We detect 12 CO emission near the optical postion of each 
QSO (offset by 1-2") and at frequencies close to those ex- 
pected for redshifted 12 CO(3-2). For each QSO, we fit a 
combination of a Gaussian and a uniform continuum to 
the spectrum corresponding to the peak signal-to-noise, 
S/N, pixel in each map and report the resulting redshift, 
line width and line flux in Table 2, as well as overplot- 
ting the fit on Fig. [2] We detect no significant continuum 
emission at 3mm in either QSO (3-<r limits given in Ta- 
ble 1), which is consistent with no contribution from syn- 
chrotron emission to the far-infrared luminosities of these 
systems. The quoted 1-a errors were obtained from a boot- 
strap error analysis on the model fit to the data. The re- 
sulting values of \ 2 indicate good fits to both emission lines: 
X 2 r = 1.0 for J0908-0034 and Xr = 1-1 for J0911±0027. Pre- 
vious results have found ~ 25% of SMGs display a double- 
peaked emission li ne, with a velocity d ifference between 
peaks ;>500 kms" 1 (|Bothwell et al.ll2012T ). indicating either 
a merger or disk-like kinematics. To test for this we attempt 
to fit a double Gaussian model to the spectra in Fig. [2] We 
find it produces a negligible improvement in x 2 , an d con- 
clude that single Gaussian fits are sufficient to describe the 
line profiles of both QSOs. 

We collapse the spectral cube over the frequency range 
corresponding to ± FWHM of the emission line to create 
the maps of the 12 CO (3-2) emission shown in Fig. [2] These 
maps are continuum subtracted, although as we noted above 
no significant continuum is detected in either source. We now 
discuss the 12 CO (3-2) emission line properties of each QSO 
in more detail. 

J 0908- 0034 — The integrated 12 CO (3-2) emission is de- 
tected at a S/N of 5.5<r. The 12 CO emission line is relatively 

4 We estimate errors by bootstrap resampling the 2SLAQ spec- 
trum, with replacement, and re-fitting a Gaussian model. 



narrow and so we bin to a frequency resolution of 10 MHz 
(~ 30 kms" 1 ), as shown in Fig. g] We find the 12 CO(3- 
2) emission line to be well fit by a single Gaussian with 
FWHM = 125 ± 25 kms" 1 , and a velocity-integrated flux 
of ScoA^ = 0.23 ± 0.05 Jy kms" 1 ; see Table 2. From the 
peak of the Gaussian distribution, we determine a CO 
redshift of z = 2.5008 ± 0.0002. We compared z C o with 
the wavelength of emission lines in the 2SLAQ spectrum of 
this QSO (Fig. [T| . We found the UV lines to be blue-shifted 
by At; ~ 100-500 km s" 1 , relative to zco, with Lya and Nv 
emission lines having the largest offset. Systematic offsets 
of UV emission lines are not unusual, and are usually at- 
tributed to resonant scattering of the e mission line, due t o 
material in the line of sight to the AGN (jSteidel et al.ll2010l ). 
At the spatial resolution of our D-configuration PdBI obser- 
vations we do not find any evidence for spatially resolved 
12 CO (3-2) emission, and channel maps show no sign of a 
velocity gradient across the source. 

J0911+0027— The integrated 12 CO emission is detected at 
a significance of 7.2<j. The 12 CO (3-2) spectrum is binned 
to a frequency resolution of 20 MHz (~ 60 km s" 1 ) , as shown 
in Fig. 11 We find the 12 CO (3-2) emission line to be well 
fitted by a single Gaussian distribution with FWHM = 
530 ± 100 kms" 1 , and a velocity integrated flux of ScoAf = 
0.62 ± 0.10 Jy kms" 1 ; see Table 2. We determine a 12 CO 
redshift of z = 2.3723 ± 0.0005 from the line fit. The rest 
frame UV emission lines in the 2SLAQ spectrum of this 
QSO are blue-shifted by Aw ~ 500-1500 kms -1 , relative to 
zco (Fig. 0. As with J0908-0034, Lya and Nv emission 
lines display the largest offsets, but strong absorption in 
Lya and Nv may affect the accuracy of the Gaussian fit to 
the emission lines. At the low spatial resolution of our D- 
configuration observations there is no evidence for spatially 
resolved 12 CO (3-2) emission, and the channel maps do not 
show any indication of velocity gradients across the line. 

We now compare our observations of these two far- 
infrared-luminous, but otherwise fairly typical, QSOs with 
the more luminous (and more massive) far-infrared bright 
QSOs previously studied, as well as with the SMG popula- 
tion. In particular, we investigate what the observed prop- 
erties of our QSOs can tell us about their relationship to 
SMGs. 

3.1 Gas Mass 

The 12 CO (3-2) line emission from each QSO pro- 
vides inform ation about t h e mo lecular gas mass of 
the system (|Solomon et al.l 1 1997T ) . We determine the 
12 C O(3-2) line luminosity, L' m /o _,i, following Eq. 3 
from ISolomon fc Vanden BoutJ (|2005 ). deriving Lc ( 3 _2) ~ 
(0.77 ±0.15) x 10 10 Kkms _1 pc 2 and Lc 0(3 -2) = (1-87 ± 
0.33) x lO^Kkms^pc 2 for J0908-0034 and J0911±0027 
respectively (Table 2). 
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Figure 3. The variation in Lco(3-2) sin 2 i/FWHM 2 , a proxy for 
molecular gas-mass fraction, as a function of far-infrared lumi- 
nosity, for high-redshift QSOs and SMGs. We adopt a mean in- 
clination for both populations, of 32.7°, appropriate for randomly 
orientated discs. For each sample we determine the median, but 
given the similarity of the results, we show only the QSO value 
(dashed line). Our proxy for gas mass fraction assumes the gas 
reservoirs in QSOs and SMGs have the same geometry, physical 
radius and mean inclination on the sky. We do not expect SMGs 
to suffer an inclination bias, but this may not be true for the 
optically-identified QSOs we study here, as the identification is 
based on their broad line properties. We indicate with an arrow 
the effect of adopting a mean inclination of 13° for the QSO pop- 
ulation, which may be more realistic, see £13.21 Our proxy also 
assumes the population have the same line brightness ratio, see 
£13,11 A value of = 0.52 is more appropriate for SMGs, and 
would further increase the separation between the implied gas 
mass fractions of SMGs and QSOs, with the later having lower 
gas fractions, and hence appearing more evolved. 

To estimate the total masses of the gas reservoirs in 
these QSOs we use the 12 CO(3-2) luminosity to estimate 
the expected 12 CO (1-0) luminosity and from that esti- 
mate the total gas mass. For th e first step we use the re- 
sults from IRiechers et al.1 (|201lh and assume a typical line 
brightness temperature ratio, for the gas in QSOs, of r3i 
= ^co (3-2) / Leo (i-o) ~ 1 0- e - thermalised gas) to con- 
vert ths_J^^O_(3-2}_lme luminosities of our QSOs, and 
the ICoppin et al.1 jjooj ) sample, to 12 CO(1-0) luminosi- 
ties. We cauti o n tha t we are extrapolating the results of 
IRiechers et al.l (|201lh to our sample of QSOs, and that in 
AGN dominated syst ems super-thermal ratios, i .e. r > 1, are 
not u ncommon (e.g. IPapadopoulos et alj|200a : llvison et al.l 
2012). We then adopt a CO-to-H2 conversion factor o f 
a = 1M ( Kkms -1 pc 2 ) _1 following IBothwell et ail (|2012h 
to convert L co (1-0) to give the gas masses (Table 2). Note 
that a is denoted as a conversion to H2, but the resulting H2 
mass is defined as the total H2+He gas mass. In this man- 
ner we derive gas masses of M gas = (0.77 ± 0.15) x 10 10 M 
and M gas = (1.9 ± 0.3) x 10 10 M for J0908-0034 and 
J091 1+0027 respectively. 

The gas masses of our far-infrared luminous QSOs 
are lower than the median gas mass of the extremely far - 
infrared luminous QSOs detected by I Coppin et al.l (|200Sl ). 
M gas = (4.3 ± 1.0) x 10 10 M Q , and t hey are also lower than 
the median gas mass of SMGs from IBothwell et al.l (|2012l ). 



Figure 4. 12 CO(3-2) line luminosity (L' co(3 _ 2) ) versus far- 
infrared luminosity (Lpijj,) for far-infrared bright QSOs and 
SMGs. We plot the best-fit line to the SMG population and we 
fix the gradient of the fit to the QSOs to that measu red for the 
SMGs, dL' co (3 _ 2) /dLpi R = 1.2 iBothwell et al.l2012t) . At a fixed 
Lfir the QSOs have L co (3—2) which is a factor of ~ 50 ± 23% 
lower than the SMGs, at a significance of 1.8-cr. If real, this differ- 
ence is in the sense that the far-infrared bright QSOs appear to 
have shorter gas consumption timescales, consistent with repre- 
senting a subsequent evolutionary phase. If we take into account 
the difference in the line brightness temperature ratios, r3i, for 
the two populations the statistical significance of this difference 
increases to 3.1a. 

M gas = (5.3 ±1.0) x 10 10 M . We note that there is signif- 
icant scatter around these averages, and the larger of our 
QSO gas masses is consistent within 3-cr. To allow a fair 
comparison, we have transformed the literature estimates of 
M gas to use a = 1. We also note that the bulk of the ob- 
servat ions come from 12 CO(3-2), although IBothwell et al.l 
(|2012l ) use r 3 i = 0.52 ± 0.09 in their conversion of L' co (3- 
2) to L co (1-0), which is appropriate for SMGs due to the 
presence of an extended reservoir of cold gas. As we will dis- 
cuss in ^4] the lower gas masses we derive for our far-infrared 
QSOs are consistent with them having recently evolved from 
a previous-SMG phase. 

3.2 Gas Fraction 

One crude measure of the evolutionary state of a system 
is the fractional contribution of molecular gas to its total 
mass. Hence, in the absence of significant replenishment of 
gas from external sources, we would expect the gas-mass 
fraction in SMGs to decline as they evolve, and thus in the 
evolutionary model we are testing, our QSOs ought to have 
lower gas-mass fractions than SMGs. 

In Fig. [3] we use a combination of observables 
which trace gas mass, Lco(3-2)> and dynamical mass, 
FWHM 2 / sin 2 i, as a proxy for gas-mass fraction, 
L CO ( 3 _2) sin 2 i/FWHM 2 . In choosing this proxy we 
make two implicit assumptions: (1) that the gas reservoir in 
both SMGs and QSOs is distributed in the same manner 
(it) the physical extent of the reservoir is similar in the 
two populations. Fig. [3] shows that the far-infrared bright 
QSOs have similar values of L CO ( 3 _ 2 ) sin 2 i/FWHM 2 on 
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average to SMGs. To determine if this means that they have 
similar gas-mass fractions (potentially indicating they are 
not more evolved) we need to consider the two simplifying 
assumptions above. 

As their sub-millimetre emission is optically thin it is 
believed that SMGs do not suffer from an inclination selec- 
tion bias and so, if we are modelling their gas distributions 
as disc- like systems jSwinbank et alj|201ll ). we can adopt a 
mean inclin ation of 32.7°, appro priate for randomly orien- 
tated discs (jBothwell et alj[2012h . Recent studies comparing 
12 CO linewidths of QSOs and SMGs found no distinction be- 
tween the populat ions, indicting that they have the same av- 
erage inclination (|Coppin et al.ll200ct ). As such, in our anal- 
ysis we adopt the same inclination for QSOs and SMGs. 
However, optically-identified QSOs such as those studied 
here may have preferentially biased orientations, as to have 
identified them as QSOs we must be able to observe their 
broad-line regions. Even considering situations where the 
obscuring torus of the AGN is not orientated with the host 
galaxy, we still expect a bias in the inclination of the gas 
reservoir. This arises from consideration of when the line of 
sight to the broad line region is directly through the disk 
of the host galaxy. If we adopt a m ean inclination of 13 °, 
which has been previously suggested (|Carilli &c Wandl2006t ). 
the QSOs sample has an ~ 85% lower gas mass fraction, on 
average, than the SMGs. The difference is significant at only 
2.0a, but is in the sense expected if QSOs are more evolved 
systems. 

To convert L^q^ 2 ) f° a total gas mass we must make 
assumptions about the temperature distribution within the 
gas. 12 CO (3-2) has an excitation temperature of 33 K and 
a critical density of 5 x 10 4 cm~ 3 , meaning it traces warm 
and/or dense environments. In QSOs it is thought that 
the gas is close to thermal equilibrium, meaning that the 
12 CO (3 -2) transition traces t he majority of t h e gas in the 
system i|Riechers et al.l l201lf ). IBothwell et al.l (|2012t ) show 
that the 12 CO(J,J-l) J>2 emission lines in SMGs are 
sub-thermally excited (see also lHarris et al.l 120101 ) , which 
is interpreted as an indication for the presence of multi- 
ple temperature components: a warm component associated 
with the star-forming regions and a cooler c omponent associ- 
ated with extended, quiescent, gas (see also lDanielson et alj 
l201ll ). Our proxy for the gas mass fraction assumes the con- 
version from Lco(3-2) t° gas mass is the same for both the 
QSO and SMG populations, which is a conservative assump- 
tion. If instead we were to include the differences in line 
brightness temperature ratio ()"3i) then the gas mass frac- 
tion is roughly doubled in SMGs, the offset between the 
populations in Fig. [3] becomes larger and the statistical for- 
mal significance of the difference increases to 3.6er (including 
a bias in the QSO inclination). We conclude that although 
the estimated gas mass fractions of far-infrared bright QSOs 
appear similar to typical SMGs, this is subject to a number 
of assumptions about the kinematics, excitation, sizes and 
orientations of the two populations. 



3.3 Gas Consumption Time Scales 

A second important measure of the evolutionary state 
of SMGs and QSOs comes from their gas-consumption 
timescales: M gas /SFR (or the inverse of this: the star- 
formation efficiency). In our evolutionary scenario the gas 
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Figure 5. A comparison of gas and SMBH masses for far-infrared 
bright QSOs and SMGs which hav e been observ e d in 12 CO (3— 
2). The QSO points are take n fromlCoppin et al, [|2008h and the 
average S MG gas mass from IBothwell et alj [|2012l h and SMBH 
mass from I Alexander et alj l|2008h . We show the predicted evolu- 
tion of the gas mass and SMBH mass in a model SMG in which 
star formation and SMBH growth occur in tandem for a period 
of 50 or lOOMyrs. Fo r the star formation w e assume a constant 
rate of SOOMQvr" 1 l|Bothwell et al]|2012l) The growth of the 
SMBH is calculated following Eq. 10 from lAlexander &: Hickoxl 
an d we assume it is accr eting with an Eddington ratio 
of rj = 0.2 (Alexander ct al. 2008). We conclude that the proper- 
ties of the far-infrared bright QSOs in our sampl e, and the QSOs 
with M B h < 1O 9 M from jCoppin et alj|2008l l. are consistent 
with their recent evolution from an earlier SMG-likc phase. The 
far-infrared bright QSOs with Mbh > 10 9 Mq are attributed as 
statistical outliers in this evolutionary scenario, as their SMBH 
masses are amongst the largest seen in the local universe. 



reservoir in SMGs is depleted by star formation as they 
transform from dusty systems to unobscured QSOs. Our 
'transition' FIR-bright QSOs are expected to represent a 
unique evolutionary phase, with prodigious star formation 
but lower gas masses, and thus shorter gas consumption 
timescales. 

The ratio of the observables, 12 CO (3-2) line luminos- 
ity, Lco(3-2)i and far-infrared luminosity, Lfir, can be used 
as a proxy for the gas consumption timescale of a galaxy. 
This ratio has the advantage that we do not hide any 
differences between populations in the assumptions about 
the conversion of 12 CO(3-2) measurement to 12 CO(1-0). 
We have therefore gathered samples of SMGs and FIR- 
bright QSOs which have been detected in 12 CO(3-2), al- 
lowing a direct comparison of observable quantities. We 
use lit erature values for the 12 CO (3-2) lumino sities of 
QSOs dCoppin et al.ll20Qgl : ISchumacher et al.ll2012l), SMGs 



(IBothwell et al 



20121) and lensed SMGs (Kneib et alj|2005l: 



Iono et al. 200E ; Swinbank et al.l 12010b ; Danielson et al.l 
201 if ) . We plot our QSOs along with these comparison sam- 
ples on Fig. [4] We see that of our observations follow a sim- 
ilar trend to the literature QSO values. Combining our re- 
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suits with the QSOs sample, we find they display ~ 50 ± 23% 
lower ratios of Lco(3-2)/Lfir than SMGs, potentially in- 
dicating shorter gas consumption timescales. The statistical 
significance of this difference is only 1.8-a, but we stress that 
this interpretation implicitly assumes the same value for r3i 
in the two populations. If we use the proposed values for 
QSOs and SMGs (see £13 - 1 p this would increase the difference 
between the two populations as the SMG gas masses will 
roughly double relative to the QSOs (significant at the 3.1- 
cr level), owing to the presence of a reservoir of colder gas in 
the SMGs, which is not traced by the 12 CO (3-2) transitio n 
and which is not seen in the QSOs (|Riechers et al.ll201lf ). 
Thus it appears that the far-infrared bright QSOs not only 
lack the extended , cool reservoir of gas seen in SMGs (e.g. 
Ilvison et al.ll20icl . l201.lh but also have less warm/dense gas 
than SMGs, suggesting that they are more evolved. 



4 DISCUSSION 

In order to study the link between QSOs and SMGs we 
have selected two far-infrared-detected QSOs at z ~ 2.5 (the 
epoch of peak activity in SMG and QSO populations) whose 
SMBH masses (~ 2 x 10 8 M©) are compar able to the aver- 
age S MBH mass for the SMG population (|Alexander et all 
l200fet ). In contrast the majority of previously 12 CO -detected 
(unlense d) far-infrared brig ht QSOs have Mbh ~ 10 9 - 
10 10 Mq (|Coppin et alj|2008t ). so much larger than the aver- 
age SMG Mbh that they cannot h ave recently evolved from 
a typical SMG (as pointed out by ICoppin et al.l l2008). It is 
precisely this mismatch in SMBH masses which our study 
addresses. Given their FIR luminosities and SMBH masses 
our new sample is a more accurate representation of the typ- 
ical far-infrared bright QSO population, and hence can be 
used to test an evolutionary link between SMGs and QSOs, 
through their gas and dynamical masses, as traced by their 
12 CO (3-2) emission. 

We show in Fig. [5] the distribution of estimated gas 
and SMBH ma sses for our two QSO s, as well as the aver- 
age SMG from IBothwell et all (I2012T) and t he QSO sample 
(including limits) from Coppin et alj (|2008T l. To investigate 
how g alaxies might evolve with time we follow ICoppin et al.l 
112003) an d overlay a simple evolutionary model based on 



gas consumption and SMBH gr owth timescales . We t ake 
the average SFR of SMGs from IBothwell et~ai1 (|2012f ) as 
5OOM yr" 1 and use this to calculate the reduction in their 
gas mass with time. In this simple model we assume a con- 
stant SFR and do not include mass loss due to winds. By 
selection our QSOs have SMBH masses similar to the aver- 
age SMGs, but to estimate the growth of the central SMBH 
on the gas consumption t imescale we follow Eq. 10 from 
lAlexander fc Hickoxl (|2012T ). In line with the measured prop- 
erties of SMGs, we limit the growth to an Eddington ratio , 
77 = 0.2 0, and an efficiency of 0.1 \ Alexander et alj|2008h . 
this then predicts the estimated growth in the SMBH mass 
in parallel to the depletion of the gas reservoir. 

In Fig. [5] we plot tracks showing the evolution in the 
expected gas and SMBH masses for the descendents of SMGs 



T\ = 0.2 indicates fast black hole growth. Within uncertainties 
it is feasible the growth is Eddington limited, jj = 1. 



after 50 and lOOMyrs. As can be seen, after lOOMyrs the 
SMGs in this model are expected to have almost completely 
depleted their massive, cold gas reservoirs and will have gas 
masses comparable to those we detect in the far-infrared 
luminous QSOs. At the same time their SMBHs will have 
grown by ~ 50%, resulting in masses similar to those seen 
in our target QSOs. As can be seen from Fig. [5j the rate of 
growth of the SMBH and the depletion of the gas reservoir 
can link the properties of a typical SMG to those seen in 
our far-infrared bright QSOs around ~ 100 Myrs later. 

Are all SMGs likely to go through a QSO-phase? If 
these three z ~ 2-3 populations (i.e. Ssso^m 5 mjy SMGs 
and QSOs with SMBH with masses of :>1O 8 M , which 
are both far-infrared-bright and -faint) are uniquely related 
through a simple evolutionary cycle, then the product of 
their respective space densities and lifetimes should be sim- 
ilar. On average we expect the SMG lifetimes to be compa- 
rable to their gas dep l etion timescales (F ig. or ~ 10 s yrs 
(|Swinbank et al.ll2006l ; iHickox et al.ll2012h a nd their volume 
densi ty at z ~ 2 - 3 is ~ 10" 5 Mpc -3 (e.g. IWardlow et al.1 
I2011D ." In comparison, the QSO volume density at z ~ 2 — 3 
with SMBH masse s above ~ 1O 8 M ) is ~ 10" 6 Mpc" 3 
Croom et ail l2009h and estimates of their lifetimes are 
>, 10 Myrs (jMartini fc Weinberdl200ll ; iHosokawal 120021 '). As 
noted by a number of previous authors, the space densi- 
ties and duty cycles of these populations are thus consis- 
tent with all SMGs subsequently transforming into a QSO. 
If they transform through a fa r-infrared bright QSO phase 
how long can t his b e? From IChapman et al.l (|2005l ) and 
IWardlow et"ail (|201ll ) we estimate that far-infrared bright 
QSOs have a volume density at z ~ 2 — 3 of ~ 10~ 7 Mpc -3 . 
By comparison to the SMG and QSO lifetimes this implies 
a duration for this transition phase of just ~ 1 Myr. Esti- 
mates of lifetimes and space densities for each population 
are uncertain at > 2x, and so a qualitative agreement is 
reasonable. If the QSO influences its host through winds 
and outflows with charac teristic velocities of ~ 1000 km s" 1 
fe.g. lrlarrison et al J2012I) then the estimated duration of the 
far-infrared bright phase would allow these winds to reach 
kpc-scales, comparable to the likely extent of the gas reser- 
voirs and star-form a tion activity withi n these systems (e.g. 
iTacconi et aLll2008l ; Ilvison et all 1201 ll ). This may help ex- 
plain the short duration of the 'transition' phase, although 
it may be exhaustion of the gas reservoir by star formation 
which is allowing these winds to propagate freely. 



5 CONCLUSIONS 

The main conclusions from our study are: 

• We have used IRAM PdBI to search for redshifted 
12 CO (3-2) emission from two far-infrared bright QSOs at 
z ~ 2.5 selected from the H-ATLAS survey. These QSOs 
were selected to have SMBH masses of Mbh < 3 x 10 s Mq , 
which are more comparable to typical SMGs than those in 
samples of high-redshift far-infrared luminous QSOs previ- 
ously detected in 12 CO. Our observations detect 12 CO (3- 
2)) emission from both QSOs and we derive line luminosi- 



ties of Lr 



L! 



CO(3-2) 



(0.77 ± 0.15) x lO^Kkms^pc 2 and 
= (1.87±0.33)xl0 10 Kkms" 1 pc 2 for J0908-0034 



CO(3-2) 



and J0911+0027 respectively (Table 2). 
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• Comparing our FIR-bright QSOs (and similar systems 
from the literature) with SMGs, we find that the QSOs have 
similar values of Lco(3-2) sm2 i/FWHM 2 , our proxy for gas 
mass fraction, to SMGs. However this is subject to a number 
of assumptions. If we consider that QSOs ha ve a biased incli- 
natio n angle of 13°, as has been suggested (|Carilli fc Wand 
120061 ). then the QSOs have an ~85% lower gas mass frac- 
tion, at a significance of 2.0-cr. In the absence of gas replen- 
ishment this is consistent with QSOs being more evolved 
systems. Furthermore, adopting the appropriate line bright- 
ness ratios for each population roughly doubles the gas mass 
in SMGs, increasing this difference further, and to a signif- 
icance of 3.6-ct. So while we see no evidence for evolution 
in the gas fraction, we require spatially-resolved 12 CO (1-0) 
observations to test this conclusively. 

• By comparing the gas consumption timescales, esti- 
mated from Lq ( 3 _2) / Lfir, for far-infrared bright QSOs 
with SMGs, we find that the QSOs have ~ 50 ± 23% lower 
consumption timescales, at a significance of 1.8-er. Adopt- 
ing appropriate line brightness ratios, of r3i = l for QSO, 
and r3i=0.52 for SMGs strengthens this conclusion to 3.1- 
cr. We conclude that far-infrared bright QSOs have a lower 
mass of warm/dense gas (probed dire ctly through 12 CO (3 - 
2)). Combined with previous results (|Riechers et al.lfeoilT ). 
showing that QSOs also lack an extended, cool reservoir 
of gas seen in SMGs, we interpret this as evidence that the 
far-infrared bright QSOs are at adifferent evolutionary stage 
than typical SMGs. In our evolutionary scenario this is con- 
sistent with far-infrared bright QSOs being in 'transition' 
from an SMG to a QSO. 

• We show that the gas and the SMBH masses in far- 
infrared bright QSOs and SMGs are consistent with a 
model where SMGs transform into QSOs on a timescale of 
~ 100 Myrs. Furthermore, the relative volume densities and 
expected durations of the SMG and QSO phases are con- 
sistent with all SMGs passing through a subsequent QSO 
phase, and we estimate that the likely duration of the far- 
infrared bright QSO phase is just ~ 1 Myr. We note that if 
necessary this duration is still sufficient to allow the QSO 
to influence the star formation and gas reservoirs across the 
full extent of the host galaxy through 1000-kms -1 winds 
and outflows. 

The scale of this study is too small to single handedly 
prove or disprove an evolutionary link between SMGs and 
QSOs. However the data we have obtained provides fur- 
ther evidence suppo rting the idea originally proposed by 
ISanders et al.l (|l988l ) that these populations are linked by 
an evolutionary sequence. We have compared 12 CO (3-2) 
detected QSOs and SMGs through a number of observable 
quantities, and find the timescales for gas depletion, and 
SMBH growth, needed to link SMGs to these sources are 
consistent. 
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